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Charmonium production and suppression in heavy-ion collisions at relativistic energies is investi- 
gated within different models, i.e. the comover absorption model, the threshold suppression model, 
the statistical coalescence model and the HSD transport approach. In HSD the charmonium disso- 
ciation cross sections with mesons are described by a simple phase-space parametrization including 
an effective coupling strength \Mi\ 2 for the charmonium states i — \c, J/ip, ip' ■ This allows to 
include the backward channels for charmonium reproduction by DD channels - which are missed 
in the comover absorption and threshold suppression model - employing detailed balance without 
introducing any new parameters. It is found that all approaches yield a reasonable description of 
J/ip suppression in S+U and Pb+Pb collisions at SPS energies. However, they differ significantly in 
the if)' /J /if) ratio versus centrality at SPS and especially at RHIC energies. These pronounced dif- 
ferences can be exploited in future measurements at RHIC to distinguish the hadronic rescattering 
scenarios from quark coalescence close to the QGP phase boundary. 

PACS numbers: 25.75.-q, 13.60.Le, 14.40.Lb, 14.65.Dw 



I. INTRODUCTION 

The dynamics of ultra-relativistic nucleus-nucleus 
collisions at Super-Proton-Synchrotron (SPS) and 
Relativistic-Heavy-Ion- Collider (RHIC) energies are of 
fundamental interest with respect to the properties of 
hadronic/partonic systems at high energy densities. Es- 
pecially the formation of a quark-gluon plasma (QGP) 
and its transition to interacting hadronic matter has mo- 
tivated a large community for more than two decades 
0,0,0,0]. However, the complexity of the dynamics has 
not been unraveled and the evidence for the formation of 
a QGP and/or the properties of the phase transition is 
much debated |5| . Apart from the light and strange flavor 
(it, u, d, d, s, s) quark physics and their hadronic bound 
states in the vacuum (p, n, tt, K, <fi, A etc.), the interest in 
hadrons with charm (c, c) has been rising continuously 
since the heavy charm quark provides an additional en- 
ergy scale, which is large compared to Aqcd- The c, c 
quark degrees of freedom are of particular interest in con- 
text with the phase transition to the QGP since cc meson 
states might no longer be formed due to color screening 

However, the suppression of J/ip and ip' mesons in 
the high density phase of nucleus-nucleus collisions at 
SPS energies 0, H E IH 111 might also be at- 
tributed to a large extent to inelastic comover scattering 
(cf. HQHElHEHIiJ andRefs. therein) pro- 
vided that the correspondin g J/^- hadron cross sections 
are of the order of a few mb El El El El El ■ The- 
oretical estimates here differ by more than an order of 
mag nitude (cf. |H El El E| 

and Refs. therein), espe- 
cially with respect to J / ip-meson scattering, such that the 
question of charmonium suppression is not yet settled. 
On the other hand, at RHIC energies further absorption 
mechanisms - such as plasma screening and gluon scat- 



tering- might play a dominant role as suggested in Refs. 
El [3(j and also lead to a substantial reduction of the 
J I if) formation in central Au+Au collisions. 

Furthermore, it has been pointed out - within statis- 
tical models - that at RHIC energies the charmonium 
formation from o pen charm and anticharm mesons might 
become essential [3lll3a . l3a . l34| and even exceed the yield 
from primary NN collisions 35] . One of the prevailing 
questions thus is if open charm mesons and charmonia 
will achieve thermal equilibrium with the light mesons 
during the nucleus-nucleus reaction. Furthermore, is the 
distribution of charm (anti)quarks over open and hid- 
den charm mesons conform with the statistical law at 
the same freeze-out parameters as anticipated in Refs. 

mEnnm. 

In fact, a previous analysis within the Hadron-String- 
Dynamics (HSD) transport model 36] has demonstrated 
that the charmonium production from open charm and 
anticharm mesons becomes essential in central Au+Au 
collisions at RHIC. This is in accordance with indepen- 
dent studies in Refs. ElE^]. On the other hand, these 
backward channels have been found to be practically neg- 
ligible at SPS energies. There is, however, an experimen- 
tal claim |37| that open charm might be enhanced by up 
to a factor of 3 in central nucleus-nucleus collisions. In 
this case the hidden charm regeneration processes might 
be essential already in (semi-)central collisions at SPS en- 
ergies [32, |33|, yj] • A possible reason of the open charm 
enhancement is an increase of the effective production 
cross sections of heavy quarks in the strongly interact- 
ing medium |3Sj . Also strong secondary meson-baryon 
channels might be responsible for this enhancement as 
pointed out in Ref. [3!j. In short, there are presently 
more open questions than solid answers. 

Here we extend our previous studies [3(| with respect 
to observable ratios of charmonium states, i.e. in partic- 
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ular the tp'/J/tp ratio which is accessable by experiment. 
We compare the HSD results to the calculations within 
the standard scenarios (including only suppression chan- 
nels) as well as within the statistical coalescence model 
(SCM) [33. 

Our work is organized as follows: In Section II we re- 
mind the reader of the standard models of charmonium 
suppression as well as of the statistical coalescence model. 
We also briefly recall the 'input' of the HSD transport 
approach with respect to charmonium and open charm 
degrees of freedom. In Section III the results of all mod- 
els are presented for S+U collisions at y/s — 20 GeV, for 
Pb+Pb collisions at y/s = 17.3 GeV and Au+Au colli- 
sions at y/s — 200 GeV. We present the yields of J/ip 
and ip' as well as their ratio as a function of centrality. 
Section IV gives a summary of our findings. 



II. DESCRIPTION OF THE MODELS 

A. The standard models 

The standard approach to charmonium production in 
heavy-ion collisions assumes that cc bound states are cre- 
ated only at the initial stage of the reaction in primary 
nucleon-nucleon collisions. During the subsequent evolu- 
tion of the system, the number of hidden charm mesons is 
suppressed by i) the absorption of pre-resonance charmo- 
nium states in nuclei (the normal nuclear suppression), 
ii) the interactions of charmonia with secondary hadrons 
(comovers) and hi) a possible dissociation of cc bound 
states in the deconfined medium. The last mechanism 
was first expected in Ref. and it was proposed, that 
charmonia might be used as a probe for deconfincmcnt 
in the state of matter created at the early stage of the 
collision. 

Two basic versions of the standard scenarios have been 
considered in the literature that both restrict to suppres- 
sion mechanisms, only. One of them, the comover model 
[l~H , assumes that the charmonium suppression increases 
gradually with the density of the strongly-interacting 
medium created in the collision. No abrupt changes of 
absorption properties of the medium take place. The 
model of Ref. [4yj represents the opposite extreme: the 
suppression sets in abruptly as soon as the energy den- 
sity exceeds a threshold value, which is a free parameter 
of this model. This version of the 'suppression-only ap- 
proach will be referred to as 'the threshold scenario'. The 
latter model is motivated by the idea that the charmo- 
nium dissociation rate is drastically higher in a quark- 
gluon-plasma (QGP) than in a hadronic medium. 

For a brief description of the 'suppression-only ap- 
proach let us consider two nuclei A and B that collide 
at impact parameter b. The number of produced hidden 
charm mesons is given by |4l( 



where af N is the production cross section of the char- 
monium species i in nucleon-nucleon (N + N) collisions, 
^A(s) is the nuclear thickness function related to the nu- 

cleon density in the nucleus, and S(b,s) < 1 is a factor 
responsible for the charmonium suppression. 

At the very initial stage charmonia experience absorp- 
tion, S = S abs , by interactions with nucleons of the col- 
liding nuclei (see, for instance, Refs. Bound cc 
states are assumed to be absorbed in the so-called 'pre- 
resonance state' before the final hidden charm mesons 
are formed. This absorption cross section is therefore 
taken to be the same for all charmonia. The cross sec- 
tion a a bs — 4.4 mb ^3 i s taken from the most recent SPS 
data analysis and is close to the theoretical prediction of 
Ref. 01 . We assume that the same cross section <r a i, s 
prevails also at RHIC energies. 

Those charmonia — that survive normal nuclear sup- 
pression — are furthermore subjected to the comover 
!QHE[I3,JllEIIi or quark-gluon plasma 
QGP) suppression [4(| . We recall that both suppression 
scenarios describe successfully the centrality dependence 
of the J/ip yield in Pb+Pb collisions at the SPS. In the 
comover approach, an additional suppression factor ap- 
pears: S = S abs S co 14], which depends on the density of 
comovers and on an effective cross section a co for charmo- 
nium dissociation by comovers. The value o~H^ = 1.0 mb 
is obtained from the fit of the NA50 data on J/tp pro- 
duction in Pb+Pb at SPS (new data were added). It 
also agrees with the NA38 data for S+U collisions. This 
value corresponds to an average over all comover species, 
relative collision energies as well as all charmonium states 
contributing to the J /ip yield through their decays. From 
a fit of the ip' to J/ip ratio in S+U collisions at SPS we 
get the value af Q — 5 mb for the effective cross section 
of ip' suppression. We will assume that the cross sections 

aio^ are the same also at RHIC energies; however, the 
charmonium suppression at RHIC becomes stronger due 
to the higher comover density. 

There are two reasons for an increased comover den- 
sity at RHIC relative to SPS: a) The multiplicity of 
produced secondary hadrons per unit rapidity interval 
at midrapidity increases by a factor of about 1.5 from 
y/s = 17 GeV to y/s = 200 GeV already in elemen- 
tary nucleon-nucleon collisions; b) the deviations from 
the wounded nucleon model become stronger at higher 
energies, which increases the comover density in central 
nucleus-nucleus collisions additionally. The centrality de- 
pendence of the number of light-flavored hadrons per unit 
pseudorapidity interval in Au+Au collisions at RHIC can 
be parametrized as [iij 



dN AuAu 



dy 



y=o 



dy 



[(l-x)N p /2 + xN coU ] 



y=o 



N AB(b) = NN 



AB / d 2 sT A (\s\)T B (\s-b\)S(b,s), (1) 



(2) 

where x = 0.11 for y/s = 200 GeV gg, N p (b) is the 
number of participants and N co u(b) is the number of col- 
lisions. Both are evaluated in the Glauber approach. 
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Calculating the centrality dependence of the charmo- 
nium suppression, it is convenient to introduce an reac- 
tant density in the plane transverse to the collision axis: 



n*(b,s)— (1 — x)n p (b, s) + 2xn c (b, s) 



(3) 



Here n p (b, s) and n c {b, s) are, respectively, the densities 
of nucleon participants and collisions in the transverse 
plane: 



N p (b) = / d 2 s n p {b,s) 



and 



N co ii{b) = d 2 s n co a(b,s). 



(4) 



(5) 



Note that the multiplicity of light-flavored hadrons is 
proportional to 



N*(b) 



d 2 s n*Jb,s). 



(G) 



Motivated by this fact, we assume that the comover den- 
sity in the transverse plane, which is needed to calculate 
S co , is proportional ton* ©. 

In contrast to the co-mover version of the suppression 
models, the threshold scenario 01 assumes that no char- 
monia are destroyed by the medium until the energy den- 
sity reaches a threshold value. The excited charmonia % c , 
which contribute about 40% to the total J ftp yield, are 
suppressed at lower energy densities compared to directly 
produced J ftp's. We have updated the fit [46| to the SPS 
data (new NA50 data 0] were added) using the cor- 
rected value of the normal nuclear absorption cross sec- 
tion a abs = 4.4±0.5 mb H^. The J ftp to Drell-Yan ratio 
in nucleon-nucleon collisions is aproximatelly the same as 



a 



NN , 



T NN 
'DY 



53. Our results are 1 



n v 



2.0 



'J/V>- 

1 1. j /ii> = 3.8 fm 2 . Here n x (nj/^) is the partici- 



inRef. |46 

fm~ 2 and nj/^, 

pant density in the transverse plane corresponding to the 
threshold energy density at which Xc-charmonia (J ftp's) 
are fully suppressed. (The change of the J ftp yield due 
to the tp' suppression is neglected.) The threshold for the 
tp' suppression n^i — 1.7 is obtained from a fit of the tp' 
to J ftp ratio in S+U collisions at SPS. 

Extrapolating to RHIC energies, one again has to take 
into account that the number of produced hadrons per 
unit rapidity and, consequently, the energy density of 
the produced medium grows with the collision energy 
and centrality. Due to deviations from the wounded 
nucleon model the charmonium suppression sets in 



1 These numbers are different from those of Ref. l48l . where an- 
other value of the J/i/i to Drell-Yan ratio in nucleon-nucleon col- 
lisions °'5y J ^/ "oy ~ 43 was assumed. 



when the effective participant density n*(b, s) © - rather 

than n p (b, s) - exceeds the threshold value. The num- 
ber of secondary hadrons per effective participant pair 
at y/s — 200 is higher than that at the SPS by a fac- 
tor of about 1.5. The critical energy density at RHIC 
is reached, therefore, at lower effective participant den- 



sity: n* = n x /1.5 ~ 1.3 fm and 



2.5 fm" 2 . 



n^/1.5. w 1.1 fm~ 2 . 



'<j/V>/1- 5 - 



Statistical coalescence model 



In contrast to the suppression scenario described in 
Section IIA, the statistical coalescence model j3l| as- 
sumes that hidden and open charm hadrons are created 
at hadronization near the point of chemical freeze-out 
which might be close to the phase boundary of the QGP. 
However, contrary to the pure thermal model |49| the to- 
tal amount of charm in the system is not assumed to be 
in chemical equilibrium. Indeed the relaxation time for 
the number of c and c is expected to exceed the lifetime 
of the system. Therefore, the total charm content of the 
final hadron system is assumed to be equal to the number 
of c and c created at the initial stage of A+A reactions 
by nucleon-nucleon collisions. Only the distribution of 
c and c among different hadron states is controlled by 
statistical laws in terms of the hadron gas (HG) model 
parameters: i.e. temperature T, baryonic chemical po- 
tential fib and volume V. It appears that the number of 
hidden charm mesons produced by a statistical coales- 
cence mechanism depends weakly on the thermodynamic 
hadronization parameters T and /ib- The charmonium 
yield is mainly defined by the average number of charmed 
quark-antiquark pairs N c£ and by the hadronization vol- 
ume parameter V. We recall that, if N cS is not much 
larger than unity, a proper account for the exact charm 
conservation becomes essential as shown in Ref. [32$ . 
This is crucial at SPS energies, where N c5 is less than 
unity, and remains essential for moderate centralities in 
Au+Au collisions at RHIC. 

The SCM formula for the total (4tt) charmonium 
yield, that takes into account exact conservation of the 
number of cc pairs, was obtained in Ref. |32| . In the 
real experimental situation, however, measurements are 
performed in a limited rapidity window Ay. In the most 
simple case, when the fraction of charmonia in the rele- 
vant rapidity window does not depend on the centrality, 
one can merely use the formula for the total yield multi- 
plied by some factor £ < 1. This approach was used in 
Refs. [33l 13 for studying the SPS data, where the mul- 
tiplicity of light hadrons, which determine the freeze-out 
volume of the system, are approximately proportional to 
the number of nucleon participants N p at all rapidities. 
At RHIC the situation is different: the total (4tt) mul- 
tiplicity of light hadrons are approximately proportional 
to the number of participants N part , while at midrapidity 
it grows faster with N part [see Eq.(|2J)]. The centrality de- 
pendence of charmonium production at different rapidi- 
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ties should, in this case, be also different. To compare 
the SCM prediction to the PHENIX data 50] , which are 
related to the J/ip yield at midrapidity dNj/^/dy, one 
has to derive a formula for the charmonium yield in a 
finite rapidity interval Ay. 

To this aim let < 1 be the probability that a c 
quark, produced in a nucleus- nucleus collision, has ra- 
pidity y within the interval Ay. The probability distri- 
bution of the number k c of c quarks inside the interval 
Ay for events with fixed total (4-7r) number N cS of cc pairs 
is given by the binomial law: 



f(k c \N c5 ) 



Nrr} 



k c \ (N c5 -k c y. ^ 



(Ay (1 - Ca«) 



JV cE -fe c 



(7) 



The probability distribution of the number k s of c's inside 
the interval Ay is assumed to be independent of k c 2 . 
It conforms to the same binomial law. Event-by-event 
fluctuations of the number of cc pairs N cS , created at the 
early stage of A + A reaction in independent nucleon- 
nucleon collisions are Poisson distributed: 



P(N c5 ; N c s) = exp (-N c5 ) 



N r .,. ! 



(8) 



The probability of cc coalescence is proportional to the 
product of their numbers and inversely proportional to 
the system volume. The proportionality coefficient de- 
pends on the thermal densities of the open and hidden 
charm hadrons, and is the same as in the case of the total 
charmonium yield 32] . 

The average multiplicity of the charmonium s pec ies i 
at fixed values of k c and k e is therefore given by |52( 



-^i ^(^c^c) ~ k c k c 



(no/2)* 



1 

VAy 



T? ■ (9) 



In deriving Eq.@ we used the fact that the thermal num- 
ber of hadrons with hidden charm is much smaller than 
that with open charm. Folding Eq.© with the binomial 
and Poisson distributions one gets 



N: 



Ay 



e Ay N c5 (n c 



1 



(no/2) 2 V A 



77- > (10) 



where no is the thermal density of all open charm 
hadrons and n\ ot is the total thermal density of the char- 
monium species i (including the decay contributions from 
the higher charmonium states). Both no and n' ot are cal- 
culated in the grand canonical ensemble with the QGP 
hadronization parameters T, //g, V Ay found from fitting 



the data of light-flavored 3 hadron yields in the rapid- 
ity interval Ay. The average number of cc pairs N c5 is, 
however, related to their total (Air) yield. 

The distinctive feature of the statistical coalescence 
model is that the ratio of multiplicities of different char- 
monium species is the same as in the equlibrium hadron 
gas. Therefore the i^ 1 to <//?/> ratio is practically indepen- 
dent on the centrality and only slightly depends on the 
collision energy (due to the change of freeze-out param- 
eters). 

In Au+Au collisions at y/s — 200 GeV the yield of 
light-flavored hadrons at midrapidity is fitted within the 
hadron gas model with T = 177 MeV and fi B = 29 MeV 
[53| . The centrality dependence of the volume is calcu- 
lated from 



V A y=l 



1 



n ch (T, hb) 



1.2 



dN 



AuAu 

ch 



drj 



(11) 



the coefficient 1.2 is needed to recalculate the number of 
particles per unit pseudorapidity (77) interval to that per 
unit rapidity (y) interval |54| . Here n c h is the charged 
hadron density calculated in the HG model. 

The average number of the initially produced cc pairs 
in our calculations is proportional to the number of bi- 
nary nucleon-nucleon collisions: 



N c5 = N coll (b)ag N /a 



NN 
inel ' 



(12) 



Our statistical coalescence model calculations are done 
under the assumption that the open charm multiplicity is 
enhanced in nucleus-nucleus collisions at the SPS accord- 
ing to the experimental claim in Ref. [371; therefore the 
effective charm production cross section a^ s N in Ea. i|12|) 
is assumed to be larger by a factor of ~ 3.5 than in ele- 
mentary nucleon-nucleon collisions. This 'enhancement' 
factor is exp ected to become weaker at larger collision 
energies [38J. Therefore, we neglect it in our calculations 
for RHIC. 

The charm production cross section, a^ z N , has been 
measured at RHIC by the PHENIX Collaboration |&5| . 
The result is consistent with PYTHIA calculations: 
a cc N ~ 650 fj,b. This gives N cE w 16—17 for central 
Au+Au collisions in line with the HSD calculations in 
Ref. 

The SCM is applicable only to large systems: N part > 
100 in Pb+Pb at the SPS EE HE 1^3- Therefore, the 
PHENIX's p + p point and the most peripheral Au+Au 
point, corresponding to N part ~ 30, cannot be used in 
the SCM fit procedure. For this reason we restrict our- 
selves only to a rough estimate of the SCM prediction for 
the charmonium yield at midrapidity at the top RHIC 
energy. 



2 This differs from Ref. 1511 . where an exact equality, k c = fcc, 
within the chosen interval Ay is assumed. In fact, the net charm 
is exactly zero only in the total system. In any finite rapidity 
interval, however, event-by-event fluctuations with k c 7^ fee are 
possible. 



3 At RHIC the strangeness as well as all other conserved charges, 
excluding charm, can be safely considered in the grand canonical 
ensemble. 
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We fix the charm production cross section in nucleon- 
nucleon collisions at its PYTHIA value, a^ 5 N — 650 /ib. 
Since there are no experimental data for the value of 
£aj/=1j one can roughly estimate it assuming approxi- 
mately the same rapidity distribution for the open charm 
and J/^'s in p + p collisions. This leads to £aj/=i ~ 0.3. 
We note that the charm rapidity distribution in Au+Au 
collisions might be broader than in p + p reactions due 
to rescattering of c and c with nucleons. This will not 
change our result essentially, however. The estimate of 
the total charm production cross section is based on the 
single electron measurement at midrapidity. Any extrap- 
olation to the total phase space has been done assum- 
ing that the charm rapidity distribution does not change 
from p+p to Au+Au. The charm production rate per bi- 
nary collision at midrapidity was found to be independent 
of the centrality (at least within the present accuracy of 
the measurement). This implies that the total charm 
production cross section should grow with the central- 
ity, if there is a broadening of the rapidity distribution. 
Both effects, the decrease of £aj/=i an d the increase of 
a^ s N nearly cancel each other in Eq. I|1(J|) such that the 
prediction of SCM does not change significantly. 



C. Open charm and charmonium dynamics in HSD 

In order to examine the dynamics of open charm 
and charmonium degrees of freedom during the forma- 
tion and expansion phase of the highly excited sys- 
tem created in a relativistic nucleus-nucleus collision 
within transport approaches, one has to know the num- 
ber of initially produced particles with core quarks, i.e. 

D, D, D*,D*,D S ,D S , D* S ,D* S , J/i>(lS), V'(2S), Xc (lP). 
In this work we follow the previous studies in Refs. 
0, 0, HiJ |5|J and fit the total charmonium production 
cross sections (i = Xc> J l^, ip') from NN collisions as a 
function of the invariant energy y/s by the function 



NN 



(s) = h 



is) (S)° fl <- 



where m, denotes the mass of charmonium i while yfsol = 
rrii + 2rriN is the threshold in vacuum. The parame- 
ters in (|13|l have been fixed to describe the J/tp and 
-0' data at lower energy (y/s < 30 GeV) as well as 
the data point from the PHENIX Collaboration jH3 
at y/s = 200 GeV, which gives a(pp ->■ J/tp + X) = 
3.99 ± 0.61(stat.) ± 0.58(sys. ± 0.40(abs)) fib for the 
total J/tp cross section. We use a = 0.2 mb, a — 
10, (3 — 0.775. The parameters fi are the fraction of 
charmonium states i. For the present study we choose 
f Xc = 0.636, fj/^p — 0.581, fyi — 0.21 in order to repro- 
duce the experimental ratio 



Xc2 



0.344±0.031 



measured in pp and ttN reactions '58, 59J as well as the 
averaged pp and pA ratio 

{B^(tP')a r )/(B^(J/tP)a J/lp ) ~ 0.0165 

(cf. the compilation of experimental data in Ref. |47jl. 
Here the experimentally measured J/tp cross section in- 
cludes the direct J/tp component (cjm) as well as the 
decays of higher charmonium states Xc,tp', i.e. 



exp 
T J/1> 



Note, we do not distinguish here the Xci(l-P) and XcsO-P) 
states. Instead, we use only the Xci(l-P) state (which we 
denote as \c), however, with an increased branching ratio 
for the decay to J/ip in order to include the contribution 
of Xc2(l-P), i.e. B( Xct -+ J/ip) = 0.54. We adopt B{^' -> 
J/tp) = 0.557 from 13- 

For the total charmonium production cross sections 
from ttN reactions we use the parametrization (in line 
with Ref. [3): 



r TTN 



{s)=fib 1 



(15) 



with 7 = 7.3 and b = 1.24 mb, which describes the ex- 
isting experimental data at low y/s reasonably well (cf. 
Fig. 3 from |5jj). y/~SQi = m,; + mjv + is the threshold 
in vacuum for ttN reactions. 

Apart from the total cross sections, we also need the 
differential distribution of the produced mesons in the 
transverse momentum px and the rapidity y (or Feyn- 
man xf) from each individual collision. We recall that 
xf = Pz/p™ ax ~ 2pz/\/~s with p z denoting the longitu- 
dinal momentum. For the differential distribution in xf 
from NN and ttN collisions we use the ansatz from the 
E672/E706 Collaboration 



dN 



dxFdpT 



(1 - \x F \) c exp(-6 PT p T ), 



(16) 



(13) where b PT = 2.08 GeV" 1 and c = a/(l + b/y/s). The 



'J I1> 



parameters a, b are choosen as a^N — 13.5, b^N = 24.9 
for NN collisions and cl^n = 4.11, b n N — 10.2 for irN 
collisions as in [36l l56| . 

The total and differential cross sections for open charm 
mesons from pp collisions, furthermore, are taken as in 
Refs. [UIH. We_thus refer to the results of Ref. H| 
which give ~16 DD pairs in central Au+Au collisions at 
yfs = 200 GeV, a factor of ^160 relative to the expected 
primordial J/tp multiplicity. 

Apart from primary hard NN collisions the open 
charm mesons or charmonia may also be generated by 
secondary meson-baryon (mB) reactions. Here we in- 
clude all secondary collisions of mesons with baryons by 
assuming that the open charm cross section (from Section 
2 of Ref. |56|') only depends on the invariant energy yfs 
and not on the explicit meson or baryon state. Further- 
more, we take into account all interactions of 'formed' 
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mesons - after a formation time of Tp — 0.8 fm/c (in 
their rest frame) |62| - with baryons or diquarks, respec- 
tively. As pointed out in Ref. the production of open 
charm pairs in central Au+Au collisions by mB reactions 
at RHIC energies is expected to be on the 10% level. 

In order to study the effect of rescattering we tenta- 
tively adopt the following dissociation cross sections of 
charmonia with baryons independent on the energy (in 
line with Refs. EUS): 

(T cSB = 6 mb; (17) 
&j/i>B = 4 mb; a XrB = 5 mb; a^ B = 5 mb. 

In (|17(l the cross section a C cB stands for a (color dipole) 
pre-resonance (cc) - baryon cross section, since the cc 
pair produced initially cannot be identified with a par- 
ticular hadron due to the uncertainty relation in energy 
and time. For the lifetime of the pre-resonance cc pair 
(in it's rest frame) a value of r cg = 0.3 fm/c is assumed 
following Ref. (63. This value corresponds to the mass 
difference of the i/j' and J/tp- 

For D, D* , D, D* - meson (tt,t], p . ui) scattering we ad- 
dress to the calculations from Ref. [22, which predict 
elastic cross sections in the range of 10-20 mb depending 
on the size of the formfactor employed. As a guideline 
we use a constant cross section of 10 mb for elastic scat- 
tering with mesons and also baryons, although the latter 
might be even higher for very low relative momenta. 

As already pointed out in the Introduction the J/ip 
formation cross sections by open charm mesons or the 
inverse comover dissociation cross sections are not well 
known and the significance of these channels is discussed 
controversely in the literature [2(| [^J |3^, 163. l65| . 
We here follow the concept of Ref. |3(| and introduce a 
simple 2-body transition model with a single parameter 
M 2 for each charmonium, that allows to implement the 
backward reactions uniquely by employing detailed bal- 
ance for each individual channel. Since the meson-meson 
dissociation and backward reactions typically occur with 
low relative momenta ('comovers') it is legitimate to 
write the cross section for the process 1 + 2 — > 3 + 4 
as 



0-1+2^3+4(5) = 2 \M, ; |Z 1 



V J Pi' 

where Ek and Sk denote the energy and spin of hadron 
k (k — 1,2,3,4), respectively. The initial and final mo- 
menta for fixed invariant energy y/s are given by 



P' 



(s — (mi + m 2 ) 2 )(s — (mi — m 2 ) 2 ) 



4s 



2 _ (g - (m 3 + m 4 ) 2 )(s - (m 3 - m 4 ) 2 ) 

Pf ~ ^ > ujj 

where m*. denotes the mass of hadron k. In (|18|) \Mi\ 2 
(i = Xo J/^iip') stands for the effective matrix element 
squared which for the different 2-body channels is taken 



of the form 

|M,| 2 = |M 4 | 2 for (tt, p) + (cc)* -► D + D (20) 
|M 4 | 2 = 3|M,| 2 for ( 7 r,p) + (cc) l ^ D* +D, 

D + D*, D* + D* 

\Mi\ 2 = ^\Mi\ 2 for {K,K*) + (cc)i^ D s + D, 

D s +D 

|M,| 2 = |M,| 2 for (K, K*) + (cc), D s + D* , 

D s + D*, D* + D, D* + D, D* + D* 

The relative factors of 3 in (|20l) are guided by the sum 
rule studies in [28| which suggest that the cross section 
is increased whenever a vector meson D* or D* appears 
in the final channel while another factor of 1/3 is in- 
troduced for each s or s quark involved. The factor 
((m 3 + rriij/^/s) 5 in (|18|l accounts for the suppression 
of binary channels with increasing yfs and has been fit- 
ted to the experimental data for the reactions tt + N — > 
p + N,iv + NJ> + N,K+ + A in Ref. H^. 

In Ref. |3jj we have used (for simplicity) the same 
matrix elements for the dissociation of all charmonium 
states i (i = Xc J/ip,ip') with mesons. However, there is 
no fundamental reason why these matrix elements should 
be identical. In the present study we will explore the 
charmonium " chemistry" explicitly and consider two dif- 
ferent scenarios - set 1: the same matrix element for all 
charmonium states i as in Ref. [3(|, and set 2: the ma- 
trix element squared for ip' is enhanced by factor of 1.5 
relative to J/4>'- 



setl 
set2 



\Mj/^\ 2 



\M Xc \ 2 



I Mo I 



(21) 



J/4,\ 



\M r \ 2 = 1.5 |M | 2 . 



We have fixed the parameter |Mo| 2 by comparison to the 
J/ip suppression data from the NA38 and NA50 Collab- 
orations for S+U and Pb+Pb collisions at 200 and 160 
A-GeV, respectively [13, 111 d (cf. Fig. d in Section 
III). We obtain the best fit for |M | 2 = 0.17 fm/GeV 2 
(which is slightly higher than in our previous study [3(| 
since the fractions of charmonium states ft have been 
also modified here). 

The advantage of the model introduced in (|18fl is that 
detailed balance for the binary reactions can be employed 
strictly for each individual channel, i.e. 



a3+4^1+2(s) = 0'l+2^3+4(s) 



(25-1 + 1) (2^ 2 
(2S 3 + l)(2S i + l) pj 



^ 4,(22) 



and the role of the backward reactions ((cc)i+meson for- 
mation by D + D flavor exchange) can be explored with- 
out introducing any additional parameter once |Mj| is 
fixed. The uncertainty in the cross sections 118|) is of 
the same order of magnitude as that in Lagrangian ap- 
proaches using e.g. SU (4) fi avor symmetry [23, H3 since 
the formfactors at the vertices are essentially unknown 
(28j. It should be pointed out that the comover dis- 
sociation channels for charmonia are described in HSD 
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with the proper individual thresholds for each channel in 
contrast to the comover absorption model described in 
Section II. A. 



We recall that (as in Refs. [H El El El El) 
the charm degrees of freedom are treated perturbatively 
and that initial hard processes (such as cc or Drell-Yan 
production from NN collisions) are 'precalculated' to 
achieve a scaling of the inclusive cross section with the 
number of projectile and target nucleons as Ap x At 
when integrating over impact parameter b. 



We typically perform 20 parallel runs for each impact 
parameter b in steps of Ab = 0.5 fm from b = 0.5 fm 
to b = 2Rt, where Rt denotes the target radius. Each 
parallel run here corresponds to a single Au+Au collision 
event. In central Au+Au collisions we have ~ 900 binary 
hard collisions at ^Jl = 17.3 GeV (cf. Fig. 8 of 0) 
and ~ 1300 at \fs = 200 GeV. In every binary collision 
we produce 1 particle for each species (i.e. J/ip, Xc i>' , 
D, D, D*,D*, D S ,D S , D*,D*), however, with a different 
weight. Thus, for 20 parallel runs we get about 1.8xl0 4 
(or 2.6xl0 4 ) perturbative particles for each species. 



For each single parallel run at fixed b we obtain the 
final particle multiplicity for all particle species as well 
as integral quantities such as the transverse energy Et 
as a function of rapidity y, the number of participants 
etc. Since we perform 20 parallel runs simultaneously, 
the spread in the distributions of particle multiplicities 
(or transverse energy) with respect to the individual runs 
provides some information on the fluctuations of particle 
multiplicities (as well as integral quantities). Vice versa, 
gating on events with fixed transverse energy Et (in an 
interval [E T - AE T /2,E T + AE T /2]) from all impact 
parameter b we obtain a distribution in the impact pa- 
rameter b that reflects the variation in centrality for the 
selected event class. However, for the observables pre- 
sented in Section III we have checked that the fluctua- 
tions in centrality have a minor impact on the normalized 
particle yields or ratios. 



The statistics is sufficiently good to reach an accuracy 
of particle yields of a few percent in central collisions. 
This accuracy becomes worth for peripheral collisions. 
Here we increase the number of parallel runs in order 
to obtain approximately the same number of charmonia 
and open charm mesons for fixed impact parameter as for 
central collisions. Note, however, that the statistics also 
becomes worth when including experimental acceptance 
cuts at SPS or RHIC energies. Thus, when comparing to 
data, the overall accuracy is only on the ± 5-7% level. 
This is also due to the fact that only some fraction of the 
initial charmonia survive the dynamical evolution due to 
a large number of dissociation reactions (see below). 
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FIG. 1: (Color online) The ratio B^o(J /ip) / 'a(DY) as a 
function of the transverse energy Et for S+U collisions at 
200 A- GeV (upper part) and Pb+Pb collisions at 160 A-GeV 
(lower part). The solid lines with the full squares indi- 
cate the HSD results, the dashed-dotted and short-dashed 
lines show the result of the suppression-only scenario (co- 
mover and threshold suppression model, respectively) while 
the long-dashed line stands for the statistical coalescence 
model (SCM). The experimental data have been taken from 
Refs. Elm 13. 



III. RESULTS FOR NUCLEUS-NUCLEUS 
COLLISIONS AT SPS AND RHIC 

A. SPS energies 

Let us compare the charmonium suppression at SPS 
energies with experimental data from the NA50 Collab- 
oration. This collaboration presents its results on J/ip 
suppression as the ratio of the dimuon decay of J/^'s 
relative to the Drell-Yan background in the 2.9 - 4.5 GeV 
invariant mass bin as a function of the transverse energy 
E Tl i.e. 



B^a(J/^)/a(DY)\ 2 . 9 . 



4.5, 



(23) 
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where B^ is the branching ratio for J /ip — > fi + [i~ . 

In the theoretical approaches we calculate the J/ip sur- 
vival probability Sj/^ defined as 



0.020 



s 



N 



J/1> 

fin 



N 



BB 



(24) 



where Ar^^ and N^g denote the final number of 
J/tp mesons and the number of J/ifi's produced ini- 
tially by BB reactions, respectively. In order to com- 
pare our calculated results to experimental data we 
need an extra input, i.e. the normalization factor 
B^^aNN{J /ip) I '&nn{DY), which defines the J/ip over 
Drell-Yan ratio for elementary nucleon-nucleon collisions. 
We choose B iui <jnn(J /tp) I<Jn n(DY) = 53, in line with a 
recent NA50 compilation |4]| obtained from experiments 
on proton collisions with lighter targets (cf. Sect. II. A.). 

The experimental ip' suppression is presented by the 
ratio 



B^J/ip -> fifi)a(J/ip)' 



(25) 



In our calculations we adopt this ratio to be 0.0165 for 
nucleon-nucleon collisions, which is again based on the 
average over pp, pd, pA reactions . 

Fig. n shows the ratio B lllli a(J/ip)/a(DY) as a func- 
tion of the transverse energy Et for S+U collisions at 200 
A-GeV (upper part) and Pb+Pb collisions at 160 A-GeV 
(lower part). The solid line gives the HSD result within 
the comover absorption scenario for the cross sections 
defined by (fTH|) while the various data points have been 
taken from Refs. [13,111113 • The dashed-dotted lines 
show results for the comover absorption scenario while 
the short dashed lines stand for the threshold suppres- 
sion model (cf. Sect. II.A.). It is seen that all models 
are compatible with the data for S+U as well as Pb+Pb, 
which is essentially due the fit of the matrix elements 
| Mi | 2 in (6) for the transport approach, a co = 1.0 nib for 
the comover model and the threshold participant densi- 
ties n x = 2.0 fm~ 2 and nj/^ — 3.8 fm -2 in the thresh- 
old suppression model. The statistical coalescence model 
(long-dashed line in the lower part of Fig. ^| also demon- 
strates a good agreement with the data for (semi) central 
(N part > 100 - -150) Pb+Pb collisions (the S+U data 
are outside its domain of applicability) due to a fit of the 
free parameters <j^ £ n and £aj,- 

For the proper description of the drop of the ratio Ij23(l 
in Pb+Pb collisions at Et ~ 100 GeV one has to take 
into account fluctuations of the transverse ene rgy 14^. l69| 
and energy losses in the dimuon event sample |3 ll |70j . To 
reproduce these effects in the transport approach, one 
would need much better statistics, which is not feasible 
at present. 

The ip'/J/ip ratio (gSJ is shown in Fig. |2 versus the 
transverse energy Et for S+U collisions at 200 A-GeV 
(upper part) and Pb+Pb collisions at 160 A-GeV (lower 
part) in comparison to the data from Refs. 0,0]. The 
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FIG. 2: (Color online) The ratio 12511 versus the transverse 
energy Et for S+U collisions at 200 A-GeV (upper part) and 
Pb+Pb collisions at 160 A- GeV (lower part). The dashed 
lines with the open circles and the solid lines with the full 
squares correspond to the HSD results calculated for two sets 
of parameters for the ip' matrix element - set 1 and set 2 12 1 1 . 
The assignment of the other lines is the same as in Fi g. U 
The experimental data have been taken from Refs. [tH l7ll|. 



dashed lines with the open circles and the solid lines with 
the full squares correspond to the HSD results calculated 
for two sets of parameters for the ip' matrix element - 
set 1 and set 2 (I2H . Here the results for 'set 1' overes- 
timate the ratio for S+U at high Et, whereas they are 
compatible with the ratio for Pb+Pb for Et > 60 GeV. 
The calculations for 'set 2' - including a larger matrix 
element for ip' - systematically lead to a lower ip'/J/ip 
ratio as a function of centrality. We note that no self 
energies for the D, D mesons have been incorporated so 
far. The latter change with baryon density and tempera- 
ture and differ for D and D mesons [f^. As pointed out 
in Ref. [73 dropping D, D masses lead to an increase of 
J/ip absorption by mesons and to a net lowering of the 
ip' I J /ip ratio for central collisions. 

The dashed-dotted lines in Fig. [3 show the results 
for the comover absorption scenario, the short dashed 



lines stand for the threshold suppression model, while the 
long dashed (constant) line indicates the SCM results for 
Pb+Pb. None of the models, however, reproduces the 
ratio for E T « 35 GcV (for Pb+Pb). All approaches 
roughly yield a constant ip'/J/ip ratio for Pb+Pb as a 
function of centrality for Et > 60 GeV. 



B. RHIC energies 

Whereas the differences between the results of the 
models are rather moderate at SPS energies due to a fit of 
the model parameters to the available data, the situation 
changes substantially at RHIC energies of y/s = 200 GeV. 
In Fig. |21 we show the calculated J/ip multiplicity per bi- 
nary collision - multiplied by the branching to dileptons 
- as a function of the number of participating nucleons 
Np ar t in comparison to the data from the PHENIX Col- 
laboration for Au + Au and pp reactions at y/s = 200 
GeV. The solid line with the full circles indicates the 
HSD results, which roughly agree with the SCM results 
(long dashed line) for N part > 100. The dashed-dotted 
line shows the results for the comover absorption scenario 
while the short dashed line stands for the threshold sup- 
pression model (with the parameters fixed at SPS). It is 
seen that the comover absorption model as well as the 
threshold model lead to an almost complete suppression 
of J I ip in central collisions (cf. also Ref. |5y|). As argued 
in Ref. [3(| this large suppression in the comover model 
is essentially due to a neglect of the backward channels 
D + D — > charmonia + meson. In fact, the HSD calcu- 
lations - that include the various backward channels - 
lead only to a moderate J/ip suppression roughly com- 
patible with the result of the SCM. This finding might 
suggest that the J/tp and open charm degrees of freedom 
reach approximate chemical equilibrium for mid-central 
and central Au+Au collisions at RHIC (see below). Un- 
fortunately, the present data from PHENIX do not allow 
to exclude any of the models so far. 

The ip 1 I J /ip ratio I125H at midrapidity provides further 
information. It is displayed in Fig. 0] for Au+Au col- 
lisions at y/s = 200 GeV versus the number of partici- 
pants N part . The dashed lines with the open circles and 
the solid lines with the full squares correspond to the 
HSD results for two sets of parameters for the ip' matrix 
element - set 1 and set 2 (|21f) . These calculations give 
the lowest ip' / J/ip ratio much below the ratio from the 
SCM (long dashed line). The geometrical comover model 
(dashed-dotted line) gives a higher ratio for central colli- 
sions than HSD. We attribute this difference to the fact 
that in the geometrical comover model only a single effec- 
tive cross section appears for all charmonia independent 
of threshold effects for individual channels. The thresh- 
old suppression model (short dashed line) provides the 
largest ip' / J/ip ratio for central collisions even above the 
SCM results. Since the predictions of the models differ 
by factors up to 4 future experiments with high statistics 
should allow to exclude at least some of them. 
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FIG. 3: (Color online) The calculated J/ip multiplicity per 
binary collision - multiplied by the branching to dileptons - 
as a function of the number of participating nucleons N par t 
in comparison to the data from the PHENIX Collaboration 
|50fl for Au + Au and pp reactions at y/s = 200 GeV. The 
assignment of the lines is the same as in Fig. flj 
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FIG. 4: (Color online) The midrapidity ratio <|25jl for Au+Au 
collisions at y/s = 200 GeV. The dashed lines with the open 
circles and the solid lines with the full squares correspond to 
the HSD results for two sets of parameters for the ip' matrix 
element - set 1 and set 2 12 U . The assignment of the lines 
is the same as in the previous figures. In addition, the line 
with crosses displays the rapidity integrated ratio 1251 from 
HSD for set 1, which is even slightly above the SCM result 
for central collisions. 



In addition, we present in Fig. ^ (by the line with 
crosses) the rapidity integrated ratio (|25|l from HSD for 
set 1, which is slightly above the SCM result for cen- 
tral collisions, however, still below the threshold model. 
This finding clearly demonstrates that midrapidity and 
rapidity integrated ratios have to be considered simul- 
taneously before conclusions on the amount of chemical 
equilibration can be drawn. 
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FIG. 5: (Color online) J/tp absorption and recreation by DD 
annihilation versus N par t for Au+Au collisions at midrapidity 
at yfs = 200 GeV. The solid line with full circles shows the 
J/tp survival probability 1241 . i.e. the number of final J/tp 
mesons over the numbers of J/tp initially produced by BB 
reactions (denoted as 'BB prod.'). Note: in all cases here the 
J/ip numbers include the J/tp from the decays of \c and tp' . 
The dashed line with full squares shows the integrated rate of 
J/tp absorption by baryons (12611 over 'BB prod.'. The dashed 
line with open triangles stands for the rate of J/tp absorption 
by mesons (1271 (normalized again to the primary 'BB prod.') 
and the dotted line with full triangles shows the integrated 
rate of J/tp recreated by DD annihilations (1281 (normalized 
to 'BB prod.'). 



C. Quantitative analysis of reactions rates from 
HSD 



Whereas the results of the HSD transport approach for 
J/tp suppression show a rough agreement with the pre- 
dictions from the statistical coalescence model for J/tp 
(cf. Fig. the tp' to J/tp ratios differ considerably at 
midrapidity (cf. Fig. f4J). This demonstrates that a full 
chemical equilibrium might not be achieved in the trans- 
port calculations since the total number of c, c quarks are 
about the same in both models. In order to understand 
these differences in more detail it is of interest to have a 
closer look at the reaction rates from the HSD approach 
in total and in a differential way with respect to rapidity. 

To this aim we show in Fig. [3] the J/tp absorption and 
recreation by DD annihilation versus the number of par- 
ticipants Np ar t for Au+Au collisions at midrapidity for 
^fs = 200 GeV. The solid line with full circles, further- 
more, shows the J/tp survival probability (|24[) . i.e. the 
number of final J/tp mesons over the numbers of J / tp ini- 
tially produced by BB reactions (denoted as 'BB prod.'). 
We mention that the J/tp numbers here include the J/tp 
from the decays of \c and ip' ■ The dashed line with full 
squares shows the integrated rate of J/tp absorption by 
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FIG. 6: (Color online) Integrated rate of J/tp (upper part), 
Xc (middle part) and tp' (lower part) absorption by mesons 
(dashed lines with open circles) in comparison to the recre- 
ation by DD annihilation (solid line with full squares) as a 
function of the impact parameter 6 for Au+Au at y^s = 200 
GeV. 



baryons 



dt- 



dN 



J/ip+B^X 

~dt 



(26) 



over the primary 'BB prod.'. As seen from Fig. [5] the 
final J/tp suppression is dominated by the dissociation 
with baryons. The dashed line with open triangles in- 
dicates the integrated rate of J / tp absorption by mesons 



dt- 



dN 



J/tp+r. 



+ D+D 



dt 



(27) 



(normalized again to the primary 'BB prod.') which 
is slightly lower than the dotted line with full triangles, 
which stands for the integrated rate of J / tp that are recre- 
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ated by DD annihilations 

J 00 dt dN ^J/^ + m m 

(normalized again to 'BB prod.'). Thus at practically 
all centralities - except for very peripheral collisions - 
the backward reactions by D + D annihilation overcom- 
pensate the 'comover' meson absorption. Nevertheless, 
both integrated rates are approximately comparable sug- 
gesting an approximate dynamical equilibrium between 
charmonia, light mesons and open charm mesons (cf. 
Ref. |_36] ) . However, a full chemical equlibrium is not 
achieved in the transport calculations since the ip 1 to J /ip 
ratio still depends on the matrix element for the charmo- 
nium+meson coupling as seen explicitly by comparing 
the results from set 1 with those from set 2. 

The question remains why the ip' / J /ip ratios at midra- 
pidity differ significantly in comparison to the SCM. To 
shed some light on this issue we show in Fig. EJthe time 
integrated rate of J/ip (upper part), Xc (middle part) and 
ip' (lower part) absorption by mesons (dashed lines with 
open circles) in comparison to the recreation by DD an- 
nihilation (solid line with full squares) as a function of 
the impact parameter b. As already seen from Fig. [5] the 
J r /ip recreation by D + D annihilation is larger than the 
J/ip dissociation with mesons. This situation is inverse 
for the \c (middle part) and ip' (lower part) and essen- 
tially related to the higher mass of Xc and ip' which lead 
to substantially larger dissociation with pions due to the 
vicinity of the D + D threshold. On the other hand the 
backward channels for Xc and ip' + meson recreation are 
suppressed by phase space relative to the channel J/ip + 
meson. One note of caution has to be added additionally: 
Due to sizeable differences in cross section the Xc and ip' 
dissociation extends to much larger times than the back- 
ward recreation channels. Thus dynamically there is no 
common freeze-out; the higher mass charmonium states 
(Xc and ip') decouple at later times than J/ip and may 
only be absorbed at late times but no longer recreated 
(see below). 

Some further information on this issue is displayed in 
Fig. {7\ where the rapidity distribution of the individual 
channels is shown for J/ip (upper part) and ip' (lower 
part). Though all production and absorption channels 
are rather flat in rapidity for -2 < y < 2 the final J/ip 
rapidity distribution (thick solid line) shows a local min- 
imum for -1 < y < 1. This effect - as a difference of 
large numbers - is related to the strong absorption with 
mesons and recreation by open charm and anticharm 
mesons. Thus the absorption of J/ip mesons relative to 
the initial production by baryon-baryon collisions (full 
squares) shows a nontrivial rapidity dependence. These 
effects are even more pronounced for the ip' (lower part) 
since the difference between the production and absorp- 
tion channels is most effectively seen around midrapidity 
-1 < 2/ < 1 , where the density of formed mesons is high 
and not very much delayed by formation time effects. 

In addition we show in Fig. [S]the ip' / J/ip ratio within 
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FIG. 7: (Color online) Rapidity distribution of J/ip (upper 
part) and ip' (calculated with set 1) production and absorp- 
tion channels for central (b= 1-4-4 fm) collisions of Au+Au 
at ,/s = 200 GeV . The ordering of the different lines is as 
follows: the solid lines with full squared stand for the rapid- 
ity distribution of J/ip (ip') mesons produced by initial BB 
collisions while the solid lines with open squared reflect the 
charmonia dissociation by baryons (B abs.); the dashed lines 
with crosses show the production by mB collisions. The dot- 
ted lines with open triangles show the J/ip(ip') dissociation 
by mesons while the dashed lines with full circles stand for 
the recreation of charmonia by D + D annihilation. The full 
solid lines give the final J/ip (upper part) and ip' (lower part) 
rapidity distributions. 



the SCM as a function of the temperature T at freeze- 
out. The midrapidity ratios from HSD for set 1 and set 
2 correspond to temperatures (in chemical equilibrium) 
of ~ 150 and 130 MeV, respectively, whereas the SCM 
default result is displayed for T=177 MeV. Thus in HSD 
the dynamical freeze-out conditions especially for ip' cor- 
respond to a later reaction phase than assumed in the 
SCM. This result is plausible in view of the large reac- 
tion cross sections for the channel D + D <-> ip' + meson. 
Note, that the rapidity integrated ratio from HSD for set 
1 corresponds to a temperature range from 175 to 190 
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FIG. 8: (Color online) The ratio if)'/J/ip within the SCM as 
a function of the temperature T. The upper and lower line of 
the shaded area show the systematic uncertainty in the ratio 
that arises from uncertainties in the branching ratios. The 
results for the midrapidity ratios from HSD for set 1 and set 
2 correspond to temperatures of ~ 150 and 130 MeV, whereas 
the SCM default result is quoted for T= 177 MeV. Note, that 
the rapidity integrated ratio from HSD for set 1 corresponds 
to a temperature range from 175 to 190 MeV. 



MeV, which is significantly higher than at midrapidity. 
Consequently, one has to consider not only midrapidity 
ratios but their rapidity dependence as well to obtain 
firm conclusions on freeze-out conditions. 



transport calculations at this energy. However, the back- 
ward D + D channels become substantial in Au + Au 
collisions at s/s = 200 GeV such that now an approxi- 
mate agreement of HSD with the statistical coalescence 
model is achieved for the J/ip suppression in mid-central 
and central collisions. We point out that a full chemical 
equilibration for the hidden and open charm degrees of 
freedom is not achieved in the transport calculations on 
the basis of hadronic interaction cross sections since the 
if)' to J/if) ratio still depends on the matrix element for 
the if)' coupling to mesons. The latter statement is solid 
since the cross sections employed for the J/ipfXc^' + 
mesons D + D channels have to be considered as upper 
limits because they are obtained from a fit to the char- 
monium suppression data from NA50 at SPS energies by 
discarding further absorption channels in a possibly pre- 
hadronic phase. 

In addition we have provided predictions for the 
tp'/J/if) ratio versus centrality, where the statistical 
coalescense model (SCM) shows a larger value than 
the HSD approach at midrapidity. On the other hand, 
rapidity integrated ratios in HSD arc slightly higher than 
the results from the SCM. This effect could be traced 
back to a significant rapidity dependence of the final ip' 
yield, since the net if)' absorption by mesons is maximal 
close to midrapidity. These pronounced differences 
can be exploited in future measurements at RHIC to 
distinguish a hadronic rescattering scenario from quark 
coalescence close to the QGP phase boundary. 



IV. SUMMARY 

In summarizing this work we have found that (in ab- 
sence of open charm enhancement in nucleus-nucleus 
collisions) the charmonium recreation by the backward 
D + D channels plays no substantial role at SPS ener- 
gies, which leads to a good agreement between the co- 
mover and threshold suppression models and the HSD 
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